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Abstract. We present a re-analysis of the simultaneous ~ 30 day lUE/XTE observation of NGC 7469 done in 
1996. Our main progress in this paper, in comparison to previous spectral analysis (Nandra et al. ^998 20001, 
is to adopt and fit directly to the data a detailed model of the Comptonized spectrum. This firstly allows to fit 
simultaneously the data from the UV to the hard X-ray band in a self-consistent way and secondly it gives direct 
constraints on the physical parameters of the disc-corona system, like the temperature and optical depth of the 
corona. 

Our results are completely consistent with a slab geometry where all the observed UV emission is supposed to 
cross the corona but more photon-starved geometries also give acceptable fits. Whatever the geometry is, the 
UV seed photon emission appears to be dominated by the reprocessing of the X-rays. 

We also found very interesting correlations between the different model parameters, the most important one being 
the anticorrelation between the corona temperature kTe and the UV flux Fuv- Such an anticorrelation is clearly 
inconsistent with a fixed disc-corona configuration and suggests a variation of the geometry of the system. We 
also find a correlation between the corona optical depth and the X-ray flux which may reflect processes linked to 
the corona formation. 

During the observations, NGC 7469 appears to accrete near its Eddington limit. This source could then be an 
example of magnetically dominated disc-corona system as recently proposed by Merloni (.2003 J . Finally, these 
data strongly support the presence of a pair free corona. 

Key words. 11.19.1; 11.09.1: NGC 7469; 13.25.2 

1. Introduction the UV bump, is produced by the reprocessing of part of 

the hot emission. Inversely the hot emission, at the origin 

The broad band opt-UV-X-7-ray spectra of Seyfert galax- ^^^^^ ^-ray component, is believed to be produced 

ies is mainly dominated by two components: an opti- ^y Compton upscattering of the soft photons, emitted by 

cal/UV bump, which was interpreted early on early as ^^jj ^^^^^^ on the coronal energetic electrons, 
the signature of cold material accreting onto a putative 

super massive black hole (Shields CnZll Malkan & Sargent The presence of a high energy cut-off near 100 

HHHl, and a broader component, covering the soft-X to keV, first detected by SIGMA and OSSE in NGC 4151 

the hard-X/soft7-ray energy band, with a rough power (Jourdain et al. lTM^ Maisack et al. lI993p . appears to be 

law shape generally cut-off above 100 keV. common in Seyfert galaxies (Zdziarski et al. Ht'^ Gondek 

These spectra are commonly interpreted in the frame- et al. 11996: Matt ■20011 Perola et al. 1200 21. The presence 

work of the reprocessing/upscattering models (e.g. Haardt of this component, and the lack of any annihilation lines 

& Maraschi. I1991L ri993|l . These models assume the pres- (while poorly constrained) are generally believed to be 

ence of two phases, a cold (generally an accretion disc) the signature of the thermal nature of the coronal plasma 

and a hot (called the corona) one, radiatively linked one (see however Malzac & Jourdain 2000), the cut-off energy 

with each other: part of the cold emission, giving birth to being then of the order of the coronal temperature. No 

cut-off but strong annihilation lines were indeed expected 

Send offprint requests to: P.O. Petrucci by the different non-thermal models develop in the 80-90s. 
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However when taking better into account the feedback 
of the accelerated particles on the acceleration process 
(Done et al.p^5(7 Henri & Pelletier lTMTl Petrucci, Henri 
& Pelletier al. i2001")l non-thermal models also predict 
spectral shape in agreement with the observations. 

Stronger and less ambiguous constraints on the nature 
of the coronal plasma are expected from variability stud- 
ies. For example, in the case of a thermal hot corona in 
radiative equilibrium with the cold phase, the X-ray spec- 
tral shape is expected to harden (the photon index F de- 
creases) when the corona temperature (and thus the high 
energy cut-off) increases (Haardt et al. I1997|l . Inversely, 
in the case of a non-thermal plasma, a softening of the 
spectrum is expected when the cut-off energy increases to 
ensure the balance between the accelerated particles and 
the shocked gas pressures (Petrucci et al. I2UU1|I . 

Such variations of the X-ray spectral shape may be 
produced by intrinsic changes of the corona properties 
(for example changes of the heating process efficiency) 
and/or by variation of the external environment like 
changes of the soft photons flux (and consequently of the 
corona cooling) produced by the cold phase. Simultaneous 
information in the UV and X-ray band are then needed 
to constrain the nature of the different processes acting 
in both phases. 

Several observations have already suggested a link 
between the UV and/or the extreme UV and the X-ray 
emission in Seyfert galaxies (Clavel et al. I19H2I Edelson 
et al. [BMI Marshall et al. ITMTI Chiang et al. IWiHI 
Uttley et al. I2()()()l Shemmer et al. I2()()l|l . showing strong 
correlations (with small or zero time lag) between the 
different fluxes, in good agreement with the reprocess- 
ing/upscattering interpretation. In some cases however, 
the X-rays and the optical/UV bands appear to be 
completely disconnected. For instance, the ~1 month 
simultaneous lUE/RXTE monitoring campaign on NGC 
7469 performed in 1996, showed only poor correlations 
between the 2-10 keV and UV fluxes (Nandra et al. 
119981 hereafter N98). Strong X-ray variations on small 
time scale were also detected but none in UV. Similar 
behavior has been observed in NGC 3516 (Edelson et al. 
I^nnni Maoz et al. HUM 120021) • These last results ruled 
out apparently, and at least for these sources, the basic 
interpretation presented above. 

The spectral analysis of the 1996 RXTE data of NGC 
7469 by Nandra et al. (2000 hereafter NOO) however 
permit to resolve these contradictions for this source. 
These authors indeed found that the X-ray spectral 
index in NGC 7469 was strongly correlated with the 
UV flux at zero time lag during the campaign. They 
interpret this result as a strong support for the X-rays 
being produced by Compton upscattering of the UV 
photons. In this framework, the observed variability of 
the UV flux is expected to directly modify the cooling 
rate of the hot corona, thus producing the observed X-ray 



spectral changes. The delays found between the RXTE 
hard X-rays with respect to the soft ones by Papadakis, 
Nandra, & Kazanas H2()()l|l also offers additional support 
to the hypothesis that the X-ray are most likely due to 
Comptonization of soft photons by hot, thermal electrons. 
The apparent disconnect between the X-ray and UV 
fluxes previously reported by N98, was then explained by 
the narrow band used in the X-ray (2-10 keV), which, 
in the presence of spectral variability, may result in 
incoherent behavior with respect to the UV one. 

The spectral treatment done by NOO assumed a 
simple power law shape for the X-ray continuum. It is 
known however that the real Comptonization spectral 
shape may be significantly different from a power law 
(Haardt HMHI Svensson [TMBI Petrucci et al. HUU^ . 
Moreover, the interpretation of the spectral variability 
in term of physical quantities, like the temperature and 
optical depth of the corona, is not straightforward with a 
phenomenological power law approximation and requires 
the use of approximate formulae. The aim of this paper 
is then to perform a re-analysis of the simultaneous 
lUE/XTE data by explicitly fitting the UV and X-ray 
continuum using a realistic Comptonization model. As 
just said, the first interest of this re-analysis is that we 
will have direct information on physical quantities char- 
acterizing the hot and cold phases. The broad band (from 
UV to hard X-rays) energy range of the Comptonization 
spectra also permit to fit simultaneously the lUE and 
RXTE data and to predict the spectral behavior at high 
energy (i.e. above 20 keV, the upper limit of our RXTE 
data) in a completely consistent way. 

The paper is organized as follows. We first present the 
data and the model used in Sect. 2. We then show, in 
Sect. 3, the results we obtain fitting simultaneously the 
lUE and XTE data and we also present the most inter- 
esting correlations we found between the different model 
parameters. We discuss these different results in Sect. 4 
and then conclude on this work in the last section. 

2. Observation and data analysis 

2.1. The data 

Many of the details of the RXTE observations are reported 
in N98 and NOO and the reader can refer to these papers for 
more details. The "L7-|-activation" model was employed 
for the background. The standard selection criteria in the 
"rex" RXTE processing script has also been applied. 

The background subtraction was however unreliable at 
high energies for sources as weak as NGC 7469 and the 
spectral analysis was therefore restricted to the 2-20 keV 
energy band. Different response matrices were calculated 
for different times of the observation using PCARSP v2.37 
but no significant differences were found between the de- 
rived spectra. We therefore used a single matrix for the 
entire observation. 
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The lUE data used here are described in Wanders et al. 
itTMTIl . NGC 7469 has also been observed by the HST FOS 
instrument in the midst of the lUE monitoring (Kriss et 
al. lHnnjl . The better quality of the HST data has allowed 
to disentangle the UV continuum more accurately from 
the broad wings of the emission lines and thus to identify 
four clean continuum windows free of any emission and/or 
absorption features, and centered on 1315, 1485, 1740 and 
1825 A. We fit the data using XSPECvll.l. Since the 
XSPEC fitting tools need spectral bin larger than 1 eV, 
we have rescaled the UV fluxes in order to have 4 windows 
of leV band width. In the following, and to be comparable 
with the results of Wanders et al. H1997|l and Kriss et al. 
H2(J(J(J|I . the UV fluxes are given in erg.s~^.cm~^.A~^ and 
are measured in a bandwidth of 21 A. 

Like in NOO, the observation have been divided into 30 
segments of approximately 1 day duration with sufficient 
photon statistics to permit precise spectral analysis. 

2.2. The model 

We fit the data using, for the continuum, the thermal 
Comptonization model in slab geometry developed by F. 
Haardt (Haardt 1993; Haardt 1994; Haardt et al.^MD- 
Hereafter we will refer to this Anisotropic Comptonization 
Code as AC2. The code treats Comptonization in steady- 
state, i.e., energy balance is implicitly assumed, so that an 
unspecified heating term balances radiative losses. Given 
the temperature and optical depth of the corona, the 
angle-dependent spectra of a disc-corona system in plane 
parallel geometry is derived using an iterative scatter- 
ing method, where the intrinsic anisotropy of the scat- 
tering is taken into account only in the first scattering 
order. The Compton recoil is also treated accurately. It 
includes also a reflection component described following 
White, Lightman & Zdziarski (,1988) and Lightman & 
White l)1988|l and assuming neutral matter. The spec- 
tral shape of the reflected photons is averaged over angles. 
It is multiplied by a first normalization factor which de- 
pends on the inclination angle (see Ghisellini, Haardt & 
Matt 119941 for details) . In addition, the usual R normal- 
ization is left free to vary in the fit procedure, so that, for 
the given inclination angle, i? = 1 corresponds to a solid 
angle, subtended by the reflector, of 27r. The resulting 
Comptonization spectral shape generally is quite different 
from the simple cut-off power law approximation generally 
used to mimic Comptonization spectra. Moreover, due to 
anisotropy effects, the first scattering order is reduced in 
comparison to the other orders, and the observed spectral 
shape is better approximated by a broken power law, the 
energy of the break depending mainly on kTg and kThh. It 
is expected to be in the 2-10 keV range for the parameter 
values we used (cf. Petrucci et al. l2()()()l for more details). 

The fit parameters in AC2 are the temperature of the 
corona kT^, its optical depth t, the temperature of the 
disc kThb (assumed as a single temperature black body) 
and the reflection normalization R. There are then two 



extra parameters compared to the simple power law -|- 
reflection model used by NOO. Other natural outputs of 
AC2 are the total (integrated over the solid angles) X-ray 
fluxes emitted by the corona (upward toward the observer 
and backward toward the disc) as well as the total UV flux 
emitted by the disc and the total flux Compton reflected 
at the disc surface. Different grids of this model are now 
available for public use as XSPEC table models^. 

It is important to point out that no links between the 
Comptonized spectrum and the soft UV disc emission are 
imposed a priori. The model simply adjusts its parameters 
to flt the data. It is only a posteriori that the resulting best 
flt values can be interpreted in a physically motivated sce- 
nario. More specifically, the code gives the ratio between 
an arbitrarily normalized input cooling flux, Fcooi, and 
the total emitted Compton spectrum. Fx- That ratio can 
then be related to a speciflc geometry of the system, e.g. 
a homogeneous plane parallel disc-corona system in radia- 
tive equilibrium as in Haardt & Maraschi H1991|l , in which 
case it must be Fcooi — Fx- This balance is imposed by 
the feed-back of the corona on the disc and viceversa in 
the absence of a dissipation mechanism intrinsic to the 
disc and is independent on the heating rate, provided it 
heats the electrons in the corona. A doubling of the heat- 
ing rate doubles Fx and Fcooi but does not change their 
ratio. If Fx> Fcooi, a photon starved conflguration is re- 
quired for instance a spherical corona. In the most general 
case, Fcooi is a fraction of the observed UV flux Fuv In 
fact, the condition Fuy — Fcooi is strictly valid only if all 
UV photons cross (and hence cool) the corona. In general 
Fuv > Fcooi- Thus, once the best fit parameters are ob- 
tained, one should seek for a self-consistent (in terms of 
physics and geometry) configuration of the disc+corona 
system. 

AC2 also does not treat the ionization and thermal 
balance of the cold phase at all. This has certainly some 
impact on the broad band spectral shape expected from 
such a model as shown by Malzac et al. (21)l)3 i , leading, for 
instance, to higher corona temperature and harder spec- 
tra in comparison to the "blackbody -I- neutral reflection" 
approximation used here. We believe however that these 
effects have no strong influence on the spectral variability 
behaviors studied in this paper. 

We have also the possibility of using a semi-spherical 
geometry for the corona. We have checked however that 
the results presented below do not signiflcantly change 
in comparison to the slab conflguration (the main change 
being the larger value of the optical depth for a given slope 
of the X-ray power law) . We thus restrict our analysis to 
the slab geometry. 

We used the uvred model for the reddening, flxing 
the extinction to the best flt value of Kriss et al. I|2000|> 
i.e. E{B — V)~ 0.12. We flxed the column density to the 
galactic one i.e. Nh — 4.8 x 10^°cm~^. For the neutral iron 
line, we simply used a gaussian (zGAUSS model of XSPEC), 
fixing the gaussian energy and width to 6.4 keV and 0.01 

^ From ,http://pitto.mib.infn.it~haardt/ATABLESj 
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Fig. 1. Data and best fit model corresponding to a soft 
(segment 2) and hard (segment 9) state of NGC 7469. The 
data are in thick (red) lines, the unfolded model (includ- 
ing the reddening and the neutral hydrogen absorption) 
in thin solid (green) lines, and the unfolded model (with- 
out reddening or neutral hydrogen absorption) in dashed 
(blue) lines. 



keV respectively. The gaussian normalization was let free 
to vary. Finally we assume a relative normalization of 1 
between the lUE and XTE data (we have checked that 
normalizations of 0.9 or 1.1 do not significantly change 
the results). 

3. Simultaneous IDE and XTE fits 

3.1. Fit results 

We have reported in Table the best fit values of the 
different model parameters, i.e. fcTg, r, fcXbf, and R, for 29 
daily segments fitting simultaneously the lUE and XTE 
data. We do not include the first segment obtained at 
the beginning of the campaign, JD-24400245.104, since 
no lUE data were available. We have also reported in the 
table the UV (at 1315 A) and X-ray (in the 2-10 keV 
range) fluxes, and the values. We have plotted in Fig. 
n]the data and the corresponding best fit models for two 
different states of NGC 7469 (a soft one corresponding to 
segment 2 and a hard one corresponding to segment 9). 

We obtained acceptable fits for each segment and the 
total = 1528 for 1334 dof. The use of the lUE and XTE 
data provide relatively good constraints on the different 
parameters of the model. Indeed, the UV data allow to 
constrain well kThb. Moreover, in Comptonization models, 
and for a given geometry, the ratio between the UV bump 
maximum and the X-ray "plateau" (i.e. the X-ray flux 
near 0.1 keV) is mainly controlled by the corona optical 
depth. Combined with the X-ray slope in the 2-10 keV 
range it then permits relatively good determination of 
r and kTf, and thus of the X-ray continuum up to the 
hard X-ray/soft 7-ray band. Consequently, the reflection 
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Fig. 2. The light curves of the different parameters of the 
model (i.e. kT^, r, kTbb and R) are plotted in the bottom 
part of the figure as well as the total (above 0.1 keV) X-ray 
flux predicted by the AC2 code (in 10~^° erg.s~^.cm~^). 
For comparison, we have also plotted, at the top, the ob- 
served UV flux at 1315A(in 10"" erg.s-^cm'^.A^^), the 
2-10 keV X-ray flux (in 10~^^ erg.s~^.cm~^) and the pho- 
ton index light curve obtained by NOO. 



component is also well constrained even if the XTE data 
permit precise spectral study only below 20 kcV. 



3.2. Light curves of fit parameters 

We have plotted the light curves of the different model 
parameters on the bottom part of Fig. |21 as well as the 
total (above 0.1 keV) X-ray flux predicted by the AC2 
code. We have also plotted the 2-10 keV X-ray and UV (at 
1315 A) light curves and the photon index one, obtained 
by NOO, on the top of the figure. 
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Table 1. Best fit values of our model parameters kTe, r, kThb and R obtained "wlien fitting the XTE and lUE data 
simultaneously. We have also reported the X-ray and UV (at 1315 A) fluxes (in erg.s^^.cm^^ and erg.s^^.cm^^.A^^ 
respectively) 



We have tested the significance of the variabihty of 
each parameter. Fitting the different light curves -with a 
constant gives a x^ for 29 dof, of 57.6, 54.9, 8.3 and 128.5 
for fcTe, T, fcTbb and R, respectively. Except fcTf,;,, the other 
parameters are thus clearly variable (significant at >99% 
confidence). 

A more sensitive way to detect parameters variations 
is the use of the F-test. For that purpose, we fit again the 
data but fixing each of the parameters (/cTg, r, fcT^b and 
R) to its weighted average value (i.e. 322 keV, 0.054, 2 eV 
and 1.1 respectively). By comparing the total of these 
fits to the when all parameters were free, we derived F- 
valucs of 4.3, 3.9, 1.2 and 4.9 for 29 additional parameters. 
Again, it impHes highly significant changes in fcTg, r and 
R (>99%) but no significant changes for fcTbf,. The black 
body temperature kT^i, is thus consistent with a constant 
and we choose to fix it, in the following, to the mean value 
of 2 eV. 

We finally remark that the total X-ray flux (above 0.1 
keV) predicted by our model also varies similarly to the 
observed 2-10 keV flux. 



3.3. Correlations 

We have reported in Table |21 the results of the linear 
(Pearson) and rank (Spearman) correlation tests between 
the different model parameters, assuming no time lags. 

The X-ray flux is apparently correlated, with signifl- 
cances >99%, with the coronal optical depth r and anti- 
correlated with R. This anticorrelation was not observed 
by NOO. We also found a correlation between the X-ray 
and line fluxes but not as strong as NOO. Concerning 
the UV flux, it appears strongly anti-correlated with the 
coronal temperature kT^. Finally, we found a strong anti- 
correlation between kTf. and r. We have plotted the most 
significant correlations in Figs. 13 and 01 i.e t and R versus 
F2-W, kTe versus -Fisis, and kT^ versus r. 

As already remarked by NOO, the cross correlation of 
light curves that have a "red noise" character (i.e. where 
the variability power scales as /~"), can produce arti- 
ficially high correlation coefficients values (Welsh 119991 
Maoz et al. l2000|l . We have thus tested the significance of 
our results by simulating a number of X-ray and UV light 
curves with "red noise" power spectra. We adopt a value 
of a of 1.3 and 1.9 for the power law index of the X-ray 
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kTe 


T 


R 


F (2-10) 


Fl315 


^ line 






-0.65 


-0.26 


-0.15 


-0.73 


-0.03 


r 


-0.79 




-0.30 


0.67 


0.37 


0.23 


R 


-0.22 


-0.21 




-0.79 


-0.05 


-0.31 


F (2-10) 


-0.20 


0.59 


-0.79 




0.35 


0.40 


Fl315 


-0.67 


0.35 


-0.07 


0.30 




-0.05 


Fline 


-0.16 


0.33 


-0.34 


0.43 


0.04 





Table 2. The Spearman (rank) correlation coefficients are reported in the upper right portion of the Table and the 
Pearson linear ones in the lower left part. The correlations have 29 points. The correlation coefficients with formal 
chance probabilities less than 1% (i.e. with correlation coefficient > 0.46) are written in bold faces. 



and UV power-density spectra respectively (cf. N98 and 
NOO). 

In 200 trials of 29 points, we never obtained Spearman 
or Pearson correlation coefficients as high as that observed 
in the real data. We thus conclude that the correlations 
shown in bold in Tableware unlikely to arise by chance. 
We discuss these different correlations in the following. 

3.3.1. kTe versus Fuv 

The coronal temperature appears strongly anticorrelated 
with the UV flux (F1315 mimics very well the total UV 
flux variations predicted by our model and in the following 
we will now only use the term Fuv instead of F1315). We 
found a rank Spearman coefficient of 0.73 while the highest 
coefficient obtained with our simulations, using 200 trials, 
was ~ 0.5 (0.6 for 400 trials). This anti-correlation clearly 
corresponds to the F vs. Fjjv correlation found by NOO. 

These authors interpret this correlation as a strong 
support for Comptonization models where any " primary" 
increase of the UV seed photons (not accompanied by an 
increase of the coronal heating rate) is generally expected 
to produce an increase of the cooling of the corona and 
thus a steepening of the spectrum and a decreasing of the 
corona temperature. We will see in the following that, in 
the context of coupled disc-corona models where the UV 
arises from reprocessing, a geometrical and/or energetic 
change of the disc-corona configuration is required. 

3.3.2. T versus F2-10 

We find a strong correlation between the coronal optical 
depth and the X-ray emission. As already noted in Sect. 
12.21 in the code we used there is no link imposed a priori 
between r and the total X-ray flux Fx and consequently 
between r and i^2-io- Such correlation is not simply ex- 
plained in the Comptonization model framework and may 
reflect a process of different nature. 

For instance, we expect that part of the X-ray flux 
produced by the corona illuminates the surrounding cold 
thick material, seed of soft photons. It is also possible that 
the hot corona is produced by e.g. evaporation of part 
of this cold matter or by magnetic buoyancy inside the 
cold matter (Meyer & Meyer-Hofmeister '1994' Meyer et 
al. 12000, Liu et al. lMJ^ Merloni & Fabian i2002i Merloni 



I2003f) . Then variations of the X-ray flux may likely modify 
the properties of the upper layers of the accretion disc, 
and thus influence the corona formation. For instance, 
thermal instabilities are known to exist in illuminated 
plasma (Rozahska & Czernv [19961 Nayakshin, Kazanas &; 
Kallman 2000), and the degree of ionization and temper- 
ature of the disc surface may rapidly increase in response 
to an increase of the illuminating flux. It may then fa- 
vor the evaporation of part of the disc matter, increasing 
consequently the corona optical depth (but see Sect. 14. 3|) . 

3.3.3. kTe vs r 

We have reported on Fig. 0] the coronal optical depth 
T versus the coronal temperature kT^ deduced from 
our fits. A strong anti-correlation is clearly visible. 
Such anti-correlation is indeed expected in the case of 
a corona-disc configuration in radiative equilibrium. It 
depends however on the real geometry of the system 
(Haardt & Maraschi. lTM^ Stern et al. ll995l Malzac et al. 
12001) . The relation expected in the case of a slab corona 
above a passive disc has been overplotted in solid line in 
Fig. ^ The data seem to be in good agreement with the 
theoretical predictions, even if they slightly stand above 
the solid line and present some dispersion. 

It is worth noting that in our fitting procedure, we 
have assumed that all the observed UV emission, what- 
ever its origin, has crossed and cooled the corona. It is 
however possible that part of the disc emission does not 
interact with the corona, like in the case of a patchy or 
spherical corona. We have simulated such a configuration 
by adding a black body component to the model used to 
fit the data. Its temperature was fixed to 2 eV and its 
normalization was fixed so that the corresponding flux (at 
1315A) was half the observed mean UV flux. In this way, 
we impose that a fixed part of the UV emission does not 
cross and cool the corona. The corresponding values of 
T and kTe have been overplotted with triangles in Fig. 
01 We omit the error bars for clarity. The triangles are 
clearly above the solid line, as expected for a more pho- 
ton starved configuration. The total that we obtained 
with this model is 1592 still for 1334 dof, slightly larger 
than without the addition of a black body component, but 
statistically equivalent. A geometry more photon-starved 
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than the simple slab corona is thus also in agreement with 
the observations. 

It is worth noting that the addition of a constant black 
body component does not substantially modify the results 
presented above. In particular, we still found correlation 
coefficients of the same order of magnitude than those 
reported in Table |21 



3.3.4. Fk^, R versus 



10 



We found only a weak correlation between the line flux and 
the X-ray flux (cf. NOO who found a stronger correlation) 
but we did find an anti-correlation, not noted by NOO, 
between R and F2-10. 

It is known that the measurement of the refiection 
component depends on the assumed model for the contin- 
uum. For example, large differences are expected between 
a simple power law shape and a real Comptonization 
model (Petrucci et al. lMIT|l . It may then explain the dif- 
ference between our results and those of NOO. It does not 
explain however why R is anticorrelated with the 2-10 keV 
fiux. Moreover these results are relatively difficult to ex- 
plain in the standard view since both the line and the 
reflection component are expected to arise from the same 
region and thus to vary in the same way. 

It is possible that part of the reflection components is 
emitted close to the central region while the other part 
may be produced by a remote reflector (indeed the pres- 
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Fig. 4. Coronal temperature vs. coronal optical depth. 
Filled circles (with errors) have been obtained assuming 
that all the observed UV flux has crossed the corona. 
Filled triangles correspond to the case where only half 
of the observed UV flux is supposed to have crossed the 
corona. The errors have been omitted in this case for clar- 
ity. The relation expected in the case of a slab corona 
above a passive disc has been overplotted in solid line (cf. 
Sect. 13.331 for details). 

ence of a distant reflector may explain the i?-i^2-io anti- 
correlation). In this case, and in the presence of flux and 
spectral variability, like for NGC 7469, the measurement 
of R and Fune may not be very reliable and meaning- 
ful (Malzac & Petrucci I2()()2|l especially with the limited 
band pass and spectral resolution of XTE and the increas- 
ing noise at high energy (near 20 keV) . In agreement with 
NOO, we thus believe that the derived values of R and Fune 
may not be truly representative of the strength and spec- 
tral shape of the real reflection components in NGC 7469. 
We will thus not discuss this parameter in the following. 
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4. Discussion 

Important results have emerged from this analysis. Firstly, 
we obtained very acceptable flts, from the UV to the X-ray 
range, of the 29 daily spectra of NGC 7469 using a real- 
istic thermal Comptonization model. This model is able 
to explain the spectral shape of the different UV/X-ray 
daily spectra but also the UV and X-ray flux variability. 
In agreement with the conclusion of NOO (see also Chiang 
20Q2 ) , the absence of a clear correlation between the UV 
and X-ray fluxes, as reported by N98, was thus clearly mis- 
leading. It simply results from the complex flux and spec- 
tral variability of this source which may produce, when 
observed in narrow energy bands, relatively unexpected 
behaviors. Berkley et al. H2000|l have failed to reproduce 
the UV and X-ray light curves of NGC 7469 using a sim- 
ple model where an X-ray point source illuminates an infl- 
nite accretion disc. This model does not consider spectral 
changes at all and thus cannot reflect the intrinsic spec- 
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tral complexity of comptonization models. We believe that 
this is the main reason of the results they obtained. 
We also note at this point that our best fit results re- 
ported in Table ^ are quantitatively different from those 
obtained by Chiang ( 2002| fitting the same data. We also 
believe that the differences are certainly due to the ap- 
proximate model used by this author to mimic thermal 
Comptonization spectra. For example, for corona optical 
depth as small as the ones he found (cf. his Fig. 5), a ther- 
mal Comptonization spectrum is formed of well separated 
bumps of different scattering orders and a simple cut-off 
power law is a really poor approximation. 

Secondly, we found very interesting correlations be- 
tween the UV/X-ray fluxes and the model parameters, 
the most important one being the ant icor relation between 
the UV flux and the corona temperature. It is the equiv- 
alent of the Fj/y-F correlation found by NOO, fitting the 
data with a simple cut-off power law model. This anti- 
correlation can be simply explained in the framework of 
thermal Comptonization models. Indeed, the increase of 
the UV soft photons flux, if " primary" , generally means 
an increase of the coronal cooling. This then produces a 
decrease of the corona temperature and a softening of the 
spectrum as observed. 

However, important results stem from this interpreta- 
tion. First, the decrease of the corona temperature due 
to the increase of the coronal cooling is completely true in 
the case of a pair-free corona. Indeed, for a pair-dominated 
corona, the increase of the cooling would correspond to 
an increase of the corona temperature rather than a de- 
crease (Ghisellini & Haardt, 11994'). As discussed in Sect. 
14.41 there are actually strong indications supporting the 
presence of a pair-free corona in NGC 7469, thus justifying 
the above interpretation. 

More importantly, we note that the Fuy-kTe anticor- 
relation cannot be reconciled with a fixed disc-corona con- 
figuration in radiative balance. As explained in Sect. 
such configuration corresponds to a fixed ratio Fx / Fuv ■ 
But there is a priori no link between Fijy and the corona 
parameters like its temperature or optical depth. fcTg and 
r have only to adjust themselves in order to insure a con- 
stant heating/cooling ratio (cf. Sect. I2.2|l . There is thus 
no reason for fcTg and/or r to change, whatever the value 
of Fuv- 

As discussed in the following sections, the Fjjv-kTe an- 
ticorrelation necessarily suggests a change of the geometry 
of the disc-corona configuration. 

4.1. UV emission dominated by reprocessing 

In the Comptonization code we used, there is no con- 
straint, a priori, on the origin of the UV emission. It 
may be dominated by the intrinsic disc emission, by 
the reprocessed one, or both may be of the same order. 
A check can be made a posteriori by comparing the 
T-kTf. values with the theoretical expectations as we 
have done in Sect. 13.3.31 and Fig. ^ If the intrinsic disc 



emission were dominant, the corona should be over- 
cooled in comparison to a corona above a passive disc, 
and the T-kT^ values should be in the lower right part 
of Fig. 01 below the solid line, which is not what we obtain. 

Another way to see this is to compare the X-ray flux 
expected to be reprocessed by the disc, Frep, and the UV 
flux which crosses and effectively cools the corona Fcooi- 
The former is simply the difference between the flux emit- 
ted by the corona backward toward the disc and the flux 
Compton reflected by the disc. For the latter, and in the 
case of a slab corona completely covering the disc, Fcooi is 
equal to the total disc emission (this does not correspond 
to the observed UV emission, since part of the disc emis- 
sion is comptonized in the corona). Both F^ep and Fcooi 
are natural outputs of the Comptonization model used to 
fit the data. We have plotted in Fig. the corresponding 
ratio Frep/ Fcooi- This ratio appears to be generally larger 
than 1. It is also variable between 1 and 1.8 and this point 
is discussed in the next section. If some part of the UV 
emission crossing the corona was intrinsic to the disc we 
would expect this ratio to be smaller than 1 which is not 
the case. 

These results thus suggest that, in the case of NGC 
7469, the intrinsic disc emission is negligible in compari- 
son to the reprocessed one. This conclusion is in complete 
agreement with that of NOO. It also agrees with the ob- 
served delays, increasing with wavelength, observed be- 
tween different UV bands and interpreted as the different 
light travel time between the illuminating X-ray source 
and the disc reprocessing regions (Wanders et al. [TM7I 
Kriss et al. 1^0001) . 

Still concerning Fig.|51 we note that, starting from the 
epoch (segm n. 15) where the ratio Frep/ Fcooi is nearly 1, 
as expected for the slab configuration, this ratio reaches 
its maxima in near correspondence with the maxima of 
the X-ray light curve (near minima for the UV fluxes). 
One can thus interpret the " oscillation" as a change in ge- 
ometry whereby the fraction of reprocessed photons which 
re-enter and cool the corona varies. 

4.2. Corona geometry variations as the origin of the 
observed variability 

Simulations of Malzac & Jourdain (|2()()l|l show that the 
time scale for a disc-corona system to reach radiative equi- 
librium, after rapid radiative perturbation in one of the 
two phases, is of the order of a few corona light crossing 
times i.e. well smaller than a day. We can thus (reason- 
ably) assume that each of the 29 daily average spectra 
and fluxes of the observation of NGC 7469 correspond to 
a disc-corona system where hot and cold phases are in 
radiative balance. 

We have also shown in the previous section that the 
disc emission is apparently dominated by the X-ray re- 
processing. In this fixed disc-corona configuration 
in radiative equilibrium should correspond to a constant 
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Fig. 5. Ratio of X-ray flux expected to be reprocessed by 
the disc, F^ep and the disc flux which crosses and actually 
cools the corona Fcooi (solid line). The theoretical predic- 
tion in the case of a slab corona above a passive disc in 
radiative equilibrium is plotted in dashed line. 



Fx/Fuv ratio (cf. Sect. 12.2(1 contrary to what we observe 
in Fig. We are thus bring to the conclusion that the 
disc-corona geometry necessarily varies during the total 
observation, producing the observed spectral and flux 
variability. 

It is not clear how the geometry really varies especially 
since we have no real constraint on the geometry itself. 
However, the fact that the X-ray emission is apparently 
the primary source of radiation (the UV being only pro- 
duced by reprocessing) , but also of variability (the total X- 
ray flux, above 0.1 keV, is indeed predicted to vary along 
the observation, cf. Fig. naturally suggests that the 
change of the disc-corona geometry is actually a change 
of the corona geometry: the larger/smaller the corona, the 
stronger/weaker the X-ray emission. 

We then propose the following scenario. Suppose that 
the origin of the day time scale variability is indeed a 
change of the global corona geometry. Suppose for exam- 
ple that the corona, which covers initially a large part 
of the disc, then becomes more patchy. Then the more 
patchy the corona, the hotter the mean corona tempera- 
ture becomes and the harder the X-ray spectrum is. At 
the same time, the corona covering factor decreasing, the 
reprocessed flux, and consequently the global UV flux, de- 
creases. This scenario thus simply explains the observed 
kTe -Fjjv anticorrelation. 

Moreover we may naturally expect the optical depth 
of the corona to depend on the disc X-ray illumination 
(which mainly occurs below the corona blobs) as already 
suggested in Sect . 13 .3!^ The variation of the X-ray flux im- 
pinging on the disc could likely modify the structure of the 
upper layers of the accretion disc and thus have an impact 
on the corona formation process. A strong X-ray emission 
would then help the evaporation of the accretion disc ma- 
terial thus increasing the corona optical depth. This could 



also explain the correlation between the corona optical 
depth and the X-ray flux. 

We of course do not assert that this scenario is unique. 
We easily imagine that the change of the disc-corona ge- 
ometry is certainly more complex, and may also imply 
changes of the disc itself. For instance, its inner region 
may become geometrically thick and optically thin, due 
to thermal instabilities (however see next section), or its 
inner radius may vary with time. Consequently, it would 
modify the UV flux entering the corona as well as the re- 
flection component. The time scale for such processes to 
occur is however expected to be quite long in comparison 
to what we observe. Anyway, whatever the scenario pro- 
posed, the present analysis strongly requires a change of 
the disc-corona configuration. On the other hand, this ap- 
pears to be a very plausible origin for the flux and spectral 
variability observed in NGC 7469. 

4.3. A magnetically dominated corona in NGC 7469? 

The mass of the black hole supposed to be present in the 
central region of NGC 7469 has been estimated in between 
10^ and 10'' solar masses. Precise measurements, using dif- 
ferent methods, favor a mass of about 7xlO^M0 (Collier 
et al. 'T9!J8' Wandel et al. I1999|l . This corresponds to an 
Eddington luminosity Ledd — 9 x 10'*'*erg.s~^. 

Our model gives us an estimate of the total luminosity 
emitted by the central engine in NGC 7469, from the UV 
to the hard X-rays (if we admit that the UV emission is 
in large part dominated by the reprocessing of X-rays, 
it is equal to the downscattered emission emitted by the 
corona). Assuming a distance of 70 Mpc (with Ho=75 
km.s~^.Mpc~^), the luminosity varies between 2 and 
3x 10^*erg.s~^ i.e. of the order of L^dd, meaning that the 
source accretes near its Eddington rate. 

It is known that for such high accretion rate, the inner 
part of the accretion disc, which is dominated by the 
radiation pressure, is viscously and thermally unstable 
(Lightman & Eardley EHZll Shakura & Sunyaev lEZHI) 
and may be subject to violent clumping instabilities 
(BegelmanEOnil Turner, Stone & Sano 2002). If the disc 
remains geometrically thin, most of the disc mass form 
very dense clumps embedded in a tenuous hot corona. In 
this case however, we expect very soft spectra unless the 
covering fraction of the cold clouds is small which results 
however in low reflection component (Malzac & Celotti 
f20D2| . This does not agree with the value of F and R 
observed in NGC 7469. 

A possible solution could be the presence of disc- 
corona configurations where a major fraction of the 
accretion power is released in the corona, as those 
proposed by Svensson & Zdziarski ifTUMjl . Merloni & 
Fabian H2002|l . and more recently by Merloni H2003|l . 
The presence of such strong coronae indeed prevents the 
development of disc instabilities. Merloni 1)2003,1 found 
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for example new thermally and viscously stable optically 
thick accretion disc solutions for magnetized turbulent 
flows. Assuming that the magnetic field amplification, via 
MRI, is balanced by buoyant escape, they obtained high 
viscosity solutions where a large fraction of the accretion 
power is liberated in the corona. These solutions appear 
only above a critical accretion rate of the order of few 
tenths of the Eddington one, which corresponds to the 
case of NGC 7469. The disc geometry being conserved in 
these solutions, we also expect the reflection normaliza- 
tion to be of the order of unity as observed. 

If this interpretation is correct, it is likely that parts 
of the disc involved in the corona formation process are 
intrinsically weakly luminous, most of the energy being 
released mechanicly/magneticaly rather than radiatively. 
The observed UV flux would then be dominated by the 
reprocessed emission in agreement with the observations. 
Moreover, it could also explain why the black body tem- 
perature kTtb deduced from our fits is significantly smaller 
than the inner temperature (at ^ 10 Schwarzschild radii) 
of ~ 50 eV expected for a standard accretion disc (i.e. 
without corona) assuming the black hole mass and accre- 
tion rate estimates of Collier et al. H1998P from optical 
data analysis. 



4.4. A pair free corona 

As already said, the code we use can give us the total 
X-ray and UV fluxes (and consequently the luminosity) 
emitted by the corona and the disc respectively. Now, from 
the variability in both bands, we can also estimate an 
upper hmit of the size of the emitting regions, Rx and 
Ruv J and consequently a lower limit on the compactnesses 
<^tLx,uv 



x,uv 



UV 



The X-rays appears to be variable on a timescale of 

less than 1 day (cf. N98, NOO) thus Ix > IT" = 77^^^ 

RinieC-^ 

where i?i is a 1 light day length. On the other hand the UV 
vary on a longer time scale of the order of 5 days (Wanders 
et al.linnZI N98; NOO; CoUier & Peterson [Mn} . Then the 



compactness ratio CR — - — > CRmin — -, — ^-^ • 

We have reported in Fig. El the (Ci?mm, ^x™) values cor- 
responding to our 29 daily observations. We obtained a 
slight correlation between these 2 variables. We also note 
that they only vary by about 50-60% during the campaign. 

These results clearly disagree with the pair-dominated 
expectations computed by Ghisellini & Haardt 1)1994(1 . 
Their Fig. 2 shows that, if the corona of NGC 7469 is 
pair dominated, its compactness ratio has to decrease and 
its high compactness to increase (both by a factor ^ 10) 
from the softest state (r=2.01, A;Te=266 keV) to the hard- 
est state (r=1.7, kTe=4:20 koV) reach by the source during 
the 29 days of observation. It would thus require very large 
changes (of factor 5 or more) of the characteristic size of 
both the UV and X-ray emitting regions for our results 
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X-ray compactness Z^*" obtained assuming a 1 light day 
for the X-ray region size (cf. Sect. 14.41 for details) 



to be compatible with the pair dominated case. This is 
in contradiction with the observations. For instance, the 
observed UV flux varying by only 50-60% imposes the size 
of the emitting region (assumed to be optically thick) to 
vary at most by the same amount. 

5. Conclusion 

We have re-analyzed in this paper the simultaneous 
30 days lUE/XTE observation of NGC 7469 done in 
1996, adopting and fitting directly to the data a detailed 
model of the Comptonized spectrum. This contrasts 
with previous spectral analysis where a simple power law 
was used to model the continuum. We were thus able 
to fit simultaneously the data from the UV to the hard 
X-ray band in a self-consistent way and to obtain direct 
constraints on the physical parameters of the disc-corona 
system, like the temperature and optical depth of the 



The model we use is able to explain well the spectral 
shape of the different UV/X-ray spectra but also the 
UV and X-ray flux variability. The absence of a clear 
correlation between the UV and X-ray fluxes, as reported 
by N98, was thus clearly misleading. It underlines the 
importance of the use of consistent models when fitting 
so intimately related energy bands. In this context, the 
observed delay of ~ 0.14 day between the UV and X-ray 
hght curves of NGC 4051, as reported by Mason et al. 
12002' may not correspond to a response delay between 
the UV and X-ray emitting region. 

Concerning NGC 7469, the data appear to be roughly 
consistent with a slab geometry where all the power is 
liberated inside the corona and where the UV flux is 
dominated by the reprocessing. We note however that 
more photon-starved geometries also agree with the data. 
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We also find very interesting correlations between 
the different model parameters as the one between the 
corona temperature and the UV flux. If it can be easily 
explained in the framework of thermal comptonization 
model (an increase of the UV flux means an increase of 
the corona cooling and thus, in general, a decrease of the 
corona temperature) , it necessarely implies slight changes 
of the disc-corona configuration during the 29 days of 
observation. We thus propose a model where the main 
origin of the variability is a change of the global corona 
geometry. The more patchy it becomes, the hotter the 
temperature. In the same time, the reprocessed flux, i.e. 
the UV emission, decreases. 

We also found an interesting correlation between the 
corona optical depth and the X-ray flux which is, a priori, 
not expected by Comptonization models and may result 
from a different, but connected, process linked to the 
corona formation. 

We show that NGC 7469 is apparently accreting at 
nearly the Eddington rate. The presence of a reflection 
component of the order of unity and of a relatively 
hard spectrum suggest the presence of a magnetically 
dominated disc-corona system in this source where most 
of the accreting energy is released in the corona. Such 
solutions have recently been shown to be thermally and 
viscously stable by Merloni H2(JU3fl . We flnally found 
strong evidences in favor of a pair free corona. 

We note finally that the lack of precise data above 
20 keV prevents however any detailed spectral study of 
the soft 7 emission of the source. For instance, the study 
of the variability of the high energy cut-off, would have 
given additionnal constraints on the model parameters, 
especially on the corona temperature. Such very broad 
band study would require simultaneous observations with 
high energy satellites hkc XMM and INTEGRAL. 
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